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Abstract Given the solution / of the sequential fractional differential equa- 
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tion aDfiaDff) + P{t)f = 0, t e [b,c], where -oo < a < b < c < +00, 
a G (|, 1) and P : [a, +00) ^ [OjPqo], -Poo < +00, is continuous, assume 
that there exist ti,t2 G [b,c] such that /(ti) = {aDtf)it2) = 0. Then, we 
estabhsh here a positive lower bound for c — a which depends solely on a, Poo- 
Such a result might be useful in discussing disco njugate fractional differential 
equations and fractional interpolation, similarly to the case of (integer order) 
ordinary differential equations. 

Key-words: Sequential fractional differential equation; Zeros of solutions; 
Fractional disco njugacy. 



1 Introduction 

The fractional calculus and its plethora of techniques for the study of frac- 
tional mathematical models in natural and social disciplines constitute al- 
ready a vast part of modern science. See the presentations from the author- 
itative monographs [21 El El El E] • The oscillation/disconjugacy theory of the 
fractional differential equations, however, is still at the very beginning, cf. 
[HI El E]- To give an example, there are no fractional counterparts of the 
fundamental comparison and separation theorems of Sturm ([TOl 17]). 

Complications such as the Faa di Bruno formula for the fractional deriva- 
tive of composite functions or the Leibniz series of the fractional derivative 
of a product P pp. 96, 98] make such enterprises extremely difficult. 

In 1918, the following striking result of W.B. Fite [2j, Theorem VII] has 
opened the way for the disconjugacy analysis [1] of (integer order) linear 
ordinary differential equations: given the real numbers a < b < c and the 
differential equation 

x" + P{t)x = 0, te[b,c\, 

where P : [a, +00) — [0, Poo], Poo < +00, is a continuous coefficient, assume 
that x{ti) = x'{t2) = for some (non-trivial) solution x. Then, 

(c — 6) ■ max{l. Poo} > 1. 

In this note, we shall establish a variant of Fife's theorem for the case 
of linear sequential fractional (ordinary) differential equations (SFDE's). In 
our result, the estimate of c — 6 is replaced with an estimate of c — a. 

To set the hypotheses, fix G (0, 1). Introduce the Riemann-Liouville 
(fractional) derivative [6, p. 68] 



t > a, 
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and the Riemann-Liouville (fractional) integral [6, p. 65] (sometimes denoted 
with Jt~^ ^ 0) 

Obviously, = / (a limit case) and = | (aA^^V)- We shall 

discuss in the next section the regularity assumptions regarding / that are 
needed in our investigation. 

Consider the following linear sequential differential equation (SFDE) 

aDtUDtf) + P{t)f = o, te[b,c], (1) 

where a E (|, l) and the coefficient P : [a, +oo) [0,Foo], -Poo < +oo, is 
assumed continuous. The equation ([T]) is equivalent to the next linear SFDE 
system 

which is a particular case of the linear SFDE system 

{aD-f)it) = Git)git) + Qit), 
UDfg){t) = R{t)f{t) + V{t), '^^''"J' 

where the coefficients G, Q, R,V : [a, +oo) —>■ M are continuous and G, R are 
bounded: \\G\\^ + \\R\\^ < +oo. 

Assume that, given the solutions {fi,gi), with « = 1,2, of the system ([3]), 
there exist the points ti G [b, c] such that 

/l(tl)=/2(tl), 9iit2)=g2it2). (4) 

Then, there will exist a positive lower bound for c — a — this is our main 
contribution in the present note. It is clear that, in the particular case of 
(12]), by taking 

(/i,^i) = (/,aA"/) and (/2,^72) = (0,0), (5) 

where / is a non-trivial solution of ([T]), we obtain the fractional version of 
Fife's theorem. 



2 Preliminaries 



To make the presentation self-contained, we shall discuss here the regularity 
of the solutions (/, g) and the transformation of ([3]) into an integral system. 
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2.1 Regularity of solutions 

Given 7 G (0, 1), consider the next integral quantity 

{Q,Mm = [^0^ds, t>a, 

where A : [a, +00) —>■ M., f : {a, +00) M are continuous and / may be 
infinite in t = a, that is 

lim(t - a)y{t) = G M. (6) 

t\a 

In the following, for such functions /, we shall use the notations 

||/IU-(a,7;M = max[(t - ay\f{t)\] 

tG|o,cJ 

and 

\\f\\L^ia,r,c)= sup[(t-an/(t)|] 

iG(a,c] 

for all b > a. 

Obviously, \\f\\L^(a,r,b,c) < ||/||L°°(a,7;c)- It is easy to see that 

\fa\ < sup ||/||L^(a,7;fe,c) = || / 1| L- (a,7;c) • (7) 
fee(a,c) 

Notice that — for a < t < c and P G (0, 1) — 

= r 1^(^)1 



^ I ii--iiL°°(a,c)||J ||L°°(a,7;c) ^^ 



= (t-a)i-'^-^S(l-7,l-/?) 

X 11^ II L-'(a,c) 11/ II L<-(a,7;c)- (8) 

We have used the change of variables s = a + X{t — a) for A G [0, 1] and the 
Beta function B, see [6j p. 6]. Remark also that 

lim(t - ay{Qp,Afm = lim(t - af{Qp,Afm = 0. 

t\a t\a 

In particular, we have the estimate 

l(Qi-.A/)(t)| < p-^^r{^ds 

< 5(1-7,7)- ||^||L-(a,c)||/||L-(a,7;c)- 
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Let us denote with C^((a, c], M) the set of functions / G C{{a, c\,M.) which 
satisfy endowed with the usual function operations. Then, it is easy to 
check that = (C^((a, c], M), || -k \\L°°{a,r,c)) is a Banach space. 

So, given / G Xy and /3 G (0, 1), the function s i — > "^(fli)^^ belongs to 
L^((a,t),M) for every t G (a, c]. 

The next estimates are crucial for establishing the Fite-type result in this 
note. 

Assume that a < b < ti < t2 < c. 

We have — via the change of variables s = a + X{t2 — a) for A G [0, 1] — 



t2 



ds 



1 {t2-sY{s-a)^ 
(t2 - af-^-^ 



a-/3-7 



dX 



1^ (1 - x)n^ 



(^2 



X 



[/1 + /2), 



a-/3-7 



(9) 



where xs denotes the characteristic function of a Lebesgue integrable subset 
S of M. For simplicity, we shall write x instead of xs from now on. 

Introduce p,q,v,w > 1 such that ^ + ^ = ^ + :^ = l and 7p, f3v < 1. 

In the following, we shall estimate the integrals Ji, I2 by means of Holder 
inequality: 



h < 



< 2 




[xix)]dx 



(10) 
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and 



h < 2^ 



< 2^ 



1 / ^ \ /3 



1-A 
dX 



x{X)dX 



(1 - A)'3- 
27 



[xwrdx 



t2-tl 



c(^;,7,/3) ■ 



t2-ti 
t-) — a 



(11) 



By taking into account (JTOl) . (flTi) and (19]), we deduce that 



(t2 - sY{s - a)T 



<(t2-a)^-^-^[c(p,/?,7) + c(t;,7,/3)] 



X max 



t2 - a J \t2- a 



= {t2-ay-^-^[cip,(3n) + civ,j,(3)]- 
Next, we would like to estimate the quantity 



where /3, 7 G (0, 1) and + 7 < 1. 
Noticing that 

5(1-7,1-/3) 
1 



t2 — ti \ ™='''(9'™> 

^2 - a 



(12) 



(tl - 0)1-/3-7 
1 



ds 



{s - ap{ti - 
*2 ds 



we get that 



ti 



1 



a {s-ay 



(ti - sY it2 - sY 



ds 



+ 



, {s-aY{ti-sY 
*^ ds 

(t2 - s)/3(s - aY 



t2 



ds 
ds 



{s-aY{t2-sY 



h + h. 
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The integral /4 was already evaluated, see ffT^ . 
We have — recall that 1 — /3 — 7>0 — 



/3 = 5(1 - 7, 1 - /?) ■ [(ti - ay-^-^ - {t, - af-^-'] < 0. 



So, 



s - an 



[{t.-sy-i {t2-sy- 



to - a 



(13) 



Finally, take / G and 7,/3 G (0, 1) such that + 7 < 1, and observe 



that 



\{Qp,Af)itl)-iQp,Af)it2)\ 

^'\Ais)-{s-arf{s)]- 



t2 



s - ap 
ds 



1 



it, - sY {t, - sf 



ds 



< \\ML°-(a,c)\\f\\L°°(a,r,b,c) 
+ 1 1 ^ I (a,c) 1 1 / 1 (a,7;c) 



ds 



{t2-sy{s-ay 



a {s-ay 



< 



(ti - sy {t2 - sY 

^^ia,c) \\f\\L^{a,r,c) " 2(^2 - ay-^~^[c{p, 7) + c{v, 7, /?)] 



ds 



X 



^2 — tl\ max{5,iD} 

t2- a 



2{t2 - a) 



1-/3-7- 



■ [c(p,/3,7) + c(t',7,/5)] 



= C-{t2~a) 



1-/3-7- 



maxjq.io} 



(14) 



where C = C(p, /3, 7) = 2[c(p, /3, 7) + c{v, 7, /?)]. 

Before going further, recall the following inequality: given the numbers 
A, 5 > and (3 G (0, 1), one has 



{A + Bf < Al^ + B^. 



(15) 
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For A = ti — a and B = t2 — ti, the inequahty ffT^ leads to 

0<{t2-af-{t,-af<{t2-t,f, 

an estimate that will be used in the next series of computations. 
Now, 

\{ti - af{Qp,Amti) - {t2 - af{Q^^Af){t2)\ 

= m-af[{Qp,Afm-{Qp,Am2)] 
-{Qp,Af)it2m2-ar-it,-an 

<{t,-af-\{Qp,Afm-{Q^,Afm\ 

+ |(QAA/)(t2)|-(t2-tlf 

< {ti - a)'^ ■ C{t2 - a) ^"^"'^^ {t2 - 

x||^||l,°°(a,c)||/||L°°(a,7;c) 

Hh - ay-P-'Bil - 7, 1 - /?) • \\A\\Loo^a,c)\\f\\L^ia,r,c) 



< 



C{c-af ■{t2-af~'^'"''^^^^^^i^ (16) 
+ (c - ay-"-^ .3(1-^,1- (3)] 
X max < (t2 — ti) ^ — ti)'^ 



< DP||i.o(,,,)||/||ioo(,,^.,) max |(t2 - ti)^^^^^, (t2 - tifj , (17) 
where = D{p, v, Pjj^b — a, c — a). We took into account the estimates dH]), 

Let us make a comment regarding D and ( fT6l) . If we don't know the sign 
of the quantity 1 — /3 — 7 ? — r then we are forced to use the raw estimate 

(^2 — a)^ ^ max{ij,™} < (fe — a)"*" ^ max{g,m} -|- (c — Ct) ^ ''^ max{g,tu} 

which shall reflect on the size of D. If, on the other hand, we have 1 — (3 — 
7 7 — r > then 

{t2 — a)"'"~^^'^^™"t9'™> < (C — a)"'"~^~'^^max{q,u,} ^ G [6, c] , 

and so 

D = C-{c- af'^~^^^^^ + (c - a)^-^-^ ■ S(l - 7, 1 - /3). (18) 

Observe that D is independent of 6, a fact of crucial importance for our 
investigation. 
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In conclusion, given / G X^, the functions t i — > {Qp,Af){t) and t i — > 
{t — a)^{Qi3^Af)(t), with < /5 < 1 — 7, are continuous in (a, c] by being 
uniformly continuous on each interval [b, c], with b G (a, c). 

In this note, by a solution of the fractional differential system ([3]) we 
mean any pair (/, g) of functions from Xi_q which verify the equations of 
the system. 

2.2 Integral representation of ([3]) 

Consider now the first of equations ([3]), namely 

{aD'^f){t) = G{t)g{t) + Q{t), te[b,c], (19) 
where a G (|, l) . 

To use the estimates from the preceding subsection, introduce the quan- 
tities 

There exist p = v > 1 such that 

-iP = f3v < ^. 

Also, q = w. 

An important consequence of this choice of constants is that 

1-/3-7 1 7 = 1 - 27 - - = - - 27 > 0, 

max|g,ty| q p 

leading to — recall (ITSj) — a quantity D which is independent of 6 G (a, c). 

27—— 27— — 

Also, since 0(^,7,7) = ^ — ^ < ^ — r- = 2^'^, we have 

(l-7p)p Q)p 

C = 2[c(p, /5, 7) + c{v, 7, /?)] = 4c(p, 7, 7) < 2'^'+^^ = 22(2-") 

and respectively 

D < 22(2-")(c-a)""^ + S(«,a)(c-a)2°~i 

< [22(2-") + 5(«,a)] i^^^ . (20) 

min{(c — a) 9 , (c — a)^~"} 



Now, recasting ( 1191) as 



LA (t-s) 



-(is 



r(i-a)[G(%(t) + g(t)], 
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we deduce that 



-ds — lim 



f{s) 



ds 



la {t-sY t\a[J^ (t-s) 

= r(l-a) f [Gis)g{s) + Qis)]ds. 

J a 

By taking into account the simple estimate 



(21) 



<B{a,l-a)- sup \{s - a)'-'' f{s) - fa\, 



the identity ( 1211) yields 



I (T^^^^ " 1 - «) = r(l - a) ^ [Gis)g{s) + Qis)]ds. 



la (t-s) 

Further, we have 



1 

r(l-a) 



^(^) M-/„S(a,l-«)- ""^^ 



a 



[G{s)g{s) + Q{s)]dsdt. 



la (x-ty- 

Next, the left hand member of the preceding identity becomes 



dt 



{x-ty~»{t- sY 



ds — faB{a, 1 — a) 



{x — ay 



a 



B{a, 1 — a) ■ / f{s)ds — faB{a, 1 — a) 



[X — a) 
a 

while the right hand member reads as, by means of the Abel computations 
P p. 32, eq. (2.13)], 

^ 'JaJa (t-s)l-« 

By differentiating with respect to x, where x > a, the identity 

(x — a)" 



B{a,l-a)- f{s)ds-faB{a,l-a) 



a 
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we get 

Notice that, since g G Xi_q,, the function s i — > G{s)g{s) + Q{s), denoted 
H, is also in Xi_a. We have Ha = G{a) ■ Qa- As /3 + 7 = (1 — a) + (1 — 
a) < 1, the integral from the right side of ([22D is {Qp^iH){x) and both 
{Qp^iH){x), {x — a)^{Qi3^iH){x) are continuous in (a,c]. 

Similarly, 

For other derivations of formulas ([221), ^^3^^ see P, p. 837 and following] 
and pp. 122-123] 



3 Main result 

Theorem 1 Set p > 1 such that (1 — a)p < | and take q = Assume 
that there exist the solutions {fi,gi) of the SFDE system ^ such that ^ 
holds for some ti, t2 ^ c]. Then, 

max {(c-a)i,(c-a)^"" I n 
m-{c-aY !^ ^> 



min{(c-a)^(c-a)l-"} 22(2-") + a) 
w/iere m = maxlHGHoo, ||-R||oo}- 
Proof. We have the relations 



f (+\ — (/')° I _J_ f* G(s)gi{£)+Q(£) , 



Take b < t < ti < c. By means of (fTTl) . we deduce that 

\{t~ay--{f,-f,m\ 

1 



[(^-a)^-"(Ql__„,G(^?2-^7l))(^) 



r(«) 

(ti-a)i-"(Qi„,,G(^2-^?i))(ti)]| 
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— -n/ N ■ II^IU°°(a,c) ||5'2 — gi\\L°°(a,l-a;c) 



X max 



D\\G\U U ( II II 

The case b < ti < t < c leads to the same conclusion. So, 

11/2 — /l||L°°(a,l-Q!;fc,c) 

^ -E' ■ IIG'lloo ■ 115^2 — 9l^L°°{a,l-a;c) 

< E ■ maxlllGlloo, \\R\\oo} ■ 1^2 - fi^llU-Ca.l-aic), (24) 

where, via 



E= max((c — a)i ,(c — aV " 
rial I 



22{2-a),^/ N max\{c-a)i,{c-a) 

<- + .{c-ar ^ , (25) 

r(") mm{{c- a)i,{c- ay-""} 

Similarly, 

115^2 gi\\L°°{a,l-a]b,c) 

^ E ■ ||-R||oo ■ 11/2 ~ fl\\L°°{a,l-a;c) 

< E ■ maxjllGlloo, ||-R||oo} ■ II/2 - /l||L°°(a,l-a;c)- (26) 

In conclusion, as E is independent of b, by taking into account I^Mj, (I26|l 
and ([7]), we get 

X<^-max{||G|U,||i?l|oc}-X, 

where X = max{||/2 - /l||L-{a,l-a;c), llfi'2 - 9l\\L^{a,l^a;c)}- 

The proof is complete. □ 



4 Particular cases 

In the particular case of ([I]),®, Theorem [T] provides an upper bound to the 
length of the open interval (a, c) in which at least one of / and aDff has no 
zero. Such an interval can be called a (fractional) disconjugacy interval of 
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SFDE ([T]). This type of investigation might be useful for problems related to 
the localization of zeros of various fractional variants of sine and cosine type 
of functions, cf. [1], P p. 19]. 

Another particular case is that of sequential relaxation-oscillation equa- 
tions, namely 

aDnaDtf)+P-f = V{t), te[b,c], 

where a G (|, l), the coefficient P > is constant and the free term V : 
[a, +oo) — > M is continuous. As before, we are interested in the solutions / 
such that /, aD^f G Xi^a- Such solutions exhibit several zeros for various 
terms V — see the numerical investigations from P, pp. 224-226, 238-242]. 
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